On the cosmic ray bound for models of extragalactic neutrino production 
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We obtain the maximum diffuse neutrino intensity predicted by hadronic photoproduction models of the type 
which have been applied to the jets of active galactic nuclei or gamma ray bursts. For this, we compare the pro- 
ton and gamma ray fluxes associated with hadronic photoproduction in extragalactic neutrino sources with the 
present experimental upper limit on cosmic ray protons and the extragalactic gamma ray background, employ- 
ing a transport calculation of energetic protons traversing cosmic photon backgrounds. We take into account the 
effects of the photon spectral shape in the sources on the photoproduction process, cosmological source evolu- 
tion, the optical depth for cosmic ray ejection, and discuss the possible effects of magnetic fields in the vicinity 
of the sources. For photohadronic neutrino sources which are optically thin to the emission of neutrons we find 
that the cosmic ray flux imposes a stronger bound than the extragalactic gamma ray background in the energy 
range between 10 5 GeV and 10 11 GeV, as previously noted by Waxman & Bahcall 1 . We also determine the 
maximum contribution from the jets of active galactic nuclei, using constraints set to their neutron opacity by 
gamma-ray observations. This present upper limit is consistent with the jets of active galactic nuclei producing 
the extragalactic gamma ray background hadronically, but we point out future observations in the GeV-to-TeV 
regime could lower this limit. We also briefly discuss the contribution of gamma ray bursts to ultra-high energy 
cosmic rays as it can be inferred from possible observations or limits on their correlated neutrino fluxes. 
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I. INTRODUCTION 

The connection between the emission of cosmic rays, 
gamma rays, and neutrinos from astrophysical accelerators is 
of considerable interest for the solution of the problem of the 
origin of cosmic rays [2], The reason why a fundamental rela- 
tion between these components must exist can be understood 
as follows. Particle acceleration mechanisms in cosmic plas- 
mas generally require the presence of a magnetic field which 
is able to confine the accelerated charged particles, i.e., elec- 
trons and protons (or ions). The accelerated electrons lose 
their energy quickly in synchrotron radiation in the magnetic 
field. These synchrotron photons provide a target for accel- 
erated protons to undergo photohadronic interactions, result- 
ing in the production of mesons, which decay. The particles 
which eventually emerge from this process are high energy 
photons, electrons (pairs), neutrons and neutrinos. Neutrinos 
are directly ejected due to their low interaction cross section. 
Gamma rays and secondary electrons initiate electromagnetic 
cascades, shifting the power from ultra-high energies to ener- 
gies below which the absorption of gamma rays by pair pro- 
duction is unimportant [3]. Finally, the neutrons, which un- 
like the protons, are not confined in the magnetic field, can 
escape and convert into cosmic ray protons after /3-decay, but 
their flux may be diminished by photoinduced n — > p conver- 
sions. The branching ratios which distribute the available en- 
ergy into the different channels are thereby generally of order 



unity. This leads to the conclusion that, cosmic proton accel- 
erators produce cosmic rays, gamma-rays and neutrinos with 
comparable luminosities [4]. 

The fundamental relation between cosmic ray and gamma 
ray production has the obvious consequence that AGN, which 
are known to produce a large fraction of the gamma-rays in the 
Universe, are a prime candidate for the sources of ultra-high 
energy cosmic rays (UHECR) [4] . The spectra of the emitted 
GeV-TeV gamma radiation of AGN also agree with the predic- 
tions of a hadronic production of these gamma rays [5]. Other 
prominent gamma-ray sources, in particular the violent events 
connected with gamma-ray bursts (GRB), have also been sug- 
gested as UHECR source candidates. Moreover, most of the 
extragalactic gamma-ray energy is found in a diffuse back- 
ground rather than in point sources, which allows for the pos- 
sibility that the UHECR sources could be relatively large ob- 
jects which would have a low gamma-ray surface brightness, 
such as radio galaxies [6], galaxy clusters [7,8], or even larger 
structures [9]. Whatever the sources are, the fundamental rela- 
tion between gamma-ray and neutrino fluxes implies that, if in 
fact the extragalactic gamma-ray emission is due to hadronic 
processes, a neutrino flux of a similar bolometric luminosity 
must exist. This prediction is the major motivation for high 
energy neutrino experiments, which are currently operated, 
under construction, or planned. In fact, most model predic- 
tions for extragalactic high-energy neutrino fluxes have been 
made by using the source model to determine the spectral 
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shape, and then by normalizing the total flux to some fraction 
of the diffuse extragalactic gamma-ray background (EGRB) 
[10,11]. 

In contrast to the limits set by gamma ray observations, the 
limits which could arise from the corresponding cosmic ray 
emission of the neutrino sources have been given little atten- 
tion. A detailed treatment of this problem regarding predic- 
tions for neutrino fluxes from the decay of topological de- 
fects (TD) has been given by Protheroe and Stanev [12], and a 
brief discussion of the possible relevance for diffuse neutrino 
fluxes from AGN by Mannheim [10]. Recently, it was pro- 
posed by Waxman & Bahcall [ 1 ] that indeed the measured flux 
of ultra-high energy cosmic rays provides the most restrictive 
limit on extragalactic diffuse neutrino fluxes for a broad class 
of sources. They claim that this cosmic ray bound is for neu- 
trinos of all energies two orders of magnitude lower than the 
bound previously used which was based on the EGRB. Be- 
sides the obvious restriction of their result to neutrinos from 
proton accelerators (i.e., excluding TD models), their claim is 
mainly based on three assumptions: (i) neutrons produced in 
photohadronic interactions can escape freely from the source, 
(ii) magnetic fields in the Universe do not affect the observed 
flux of extragalactic cosmic rays, and (iii) the overall injec- 
tion spectrum of extragalactic cosmic rays is ctE~ 2 . A key 
role is played by their assumption (iii): By assuming a specific 
cosmic ray input spectrum, they can normalize their bound 
at the ultra high energies, where they argue that also (ii) ap- 
plies. Assumption (i) is justified by showing that some partic- 
ular sources of specific interest, like the TeV-blazar Mrk501, 
or also GRB, are transparent to the emission of neutrons. The 
authors claim that this new bound set by cosmic ray data es- 
sentially rules out the hypothesis that hadronic processes in 
AGN jets can produce the EGRB, and consequently that their 
neutrino fluxes are overestimated. 

The purpose of the present paper is to revisit the role cos- 
mic ray observations can play to constrain models of neutrino 
production. The paper is organized as follows. In Section II, 
we give a brief review on the properties of photohadronic in- 
teractions, and derive the production spectra for cosmic rays 
and neutrinos for power-law photon target spectra. In Sec- 
tion III we briefly describe the effect of the propagation of 
extragalactic cosmic rays. We then follow Waxman & Bah- 
call in deriving a cosmic ray bound on neutrino fluxes, adopt- 
ing their assumptions (i) and (ii), but instead of assuming a 
specific cosmic ray injection spectrum, we assume a specific 
spectrum for the observable extragalactic cosmic ray flux, 
which is constructed such that it complies with all existing ob- 
servational limits on the cosmic ray proton intensity. In Sec- 
tion IV we turn to AGN, and discuss in particular the photo- 
hadronic opacity of blazar jets as can be estimated from obser- 
vations. We shall show that most AGN in fact have large pho- 
tohadronic opacities at ultra-high energies, and we derive an 
upper bound for the neutrino contribution from AGN consid- 
ering the effect of neutron opacity. In Section V, we discuss 
the possible effect on cosmic ray protons from neutron decay 
of the magnetic fields known to exist in clusters of galaxies, 
and radio galaxies which are considered the hosts of gamma- 
ray emitting AGN. We derive critical energies below which 



they could increase the bound. We conclude by discussing the 
combined effect of our results, and to what extent the cosmic 
ray data can indeed constrain models for expected neutrino 
fluxes, and vice versa. 



II. COSMIC RAY, GAMMA-RAY, AND NEUTRINO 
EMISSION FROM EXTRAGALACTIC PROTON 
ACCELERATORS 

In this section, we obtain the form of the spectra of cosmic 
rays, gamma-rays, and neutrinos escaping from cosmic pro- 
ton accelerators. Here, we shall assume that protons are con- 
fined within the acceleration region whereas photoproduced 
neutrons may escape to become cosmic rays. This is justified 
in particular for cosmic ray sources connected to relativistic 
outflows, like AGN jets or GRB, since the life time of pro- 
tons is here limited by adiabatic energy losses and generally is 
much shorter than the diffusive escape time. For example, if 
protons are accelerated near the beginning of an AGN jet, say 
where the jet width is ^10 16 cm and where the gamma-rays 
are probably produced, and are released near the end of the 
jet after it's width has expanded to at least a few parsec, then 
their energies on release will be down by a least two orders 
of magnitude. As a consequence, only protons resulting from 
the decay of neutrons which have escaped from the acceler- 
ation region contribute significantly to the cosmic ray spec- 
trum. Escaping neutrons are produced by accelerated protons 
which interact with ambient soft photons, together with pions 
which decay into neutrinos. 

We shall first discuss the properties of photohadronic in- 
teractions. We shall then obtain the ambient proton spectrum 
in the emission region resulting from shock acceleration fol- 
lowed by radiative cooling, adiabatic losses, and advection 
away from the shocked region where the photon density is 
sufficient for efficient photoproduction of neutrons. From this 
we obtain the form of the spectra of neutrons and neutrinos on 
production, and the escaping neutron spectrum which may be 
modified by neutron absorption in photohadronic n — > p con- 
versions. 

A. Photohadronic interactions 

Photohadronic interactions can be divided into two pro- 
cesses: photoproduction of pions (and other mesons), and 
Bethe-Heitler production of e ± pairs. Charged pions decay 
as 7r ± — > ^Vfj,, ^ ± — > e ± v fl v e (here and in the follow- 
ing we disregard the difference between neutrinos and anti- 
neutrinos), neutral pions decay into gamma-rays as n° — ► 77. 
Electrons and positrons (from pion decay and Bethe-Heitler 
production) cascade in the magnetic field and radiation field, 
and so can be assumed to convert all their energy into syn- 
chrotron radiation in the magnetic field required for the accel- 
eration of protons. The production of charged pions allows the 
production of secondary neutrons through isospin exchange. 

The physics of photohadronic interactions in ambient pho- 
ton spectra has been extensively studied in Monte-Carlo simu- 
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lations [13,14]. The properties of the production of secondary 
particles can be expressed in the fraction £ of the proton en- 
ergy given to a specific particle component per interaction. 
For neutrinos, gamma-rays, and neutrons, the values 

0.1 (1) 

£n « 0.5 , (2) 

respectively, have been found for power law target spectra typ- 
ical in AGN jets, while for GRB target spectra the values are 
£„ w £ 7 « £„ « 0.2 [15]. The energy per particle in units of 
the proton energy have been found for neutrinos and neutrons 
as (E v ) /E p « 0.033 and (£„} /S p w 0.83, respectively, 
while for GRB they are {E v ) / E p ~ 0.02 and (E n ) / E p ~ 
0.5, respectively [15]. From this we can immediately define 
the relative energy of escaping neutrinos and neutrons as 

fbm = (E v ) I (E n ) « 0.04 (3) 

for both AGN and GRB target spectra. The fractional energy 
loss of the proton per interaction is k p w 0.2 in the AGN 
case, and k p m 0.5 for GRB [15]. Note that the quantities 
given here as typical for GRB apply only at ultra-high proton 
energies, at lower energies they approach the values found for 
AGN [15]. 

Electromagnetic radiation initiated by the Bethe-Heitler 
pair production, and by photons and electrons from neutral 
and charged pion decay, are reprocessed in synchrotron pair 
cascades. This energy will emerge as a component of the 
gamma ray background radiation, for AGN mainly in an en- 
ergy range lOMeV— 1 TeV. The contribution of the Bethe- 
Heitler process to the production of gamma-rays depends on 
the target energy spectrum index, a, since its cross section 
peaks at energies about two orders of magnitude lower than 
that of photopion production. Assuming that the power law 
extends over this range without change, one can find the rela- 
tion 

L 7 = [1 + cxp(5a - 5)]L V , (4) 

where L 7 and L v are the bolometric photohadronic luminosi- 
ties in gamma rays and neutrinos [16], and we note that for 
a = 1, L 1 = 2L V . We also note that in general, L 7 > 
L v holds as a direct consequence of the isospin-symmetry of 
charged and neutral pions — hence, for any kind of neutrino 
production involving pion decay, the bolometric flux in corre- 
lated photons sets an robust upper limit on the possible bolo- 
metric neutrino flux. 



B. Ambient proton spectrum 

We assume a spectrum of protons on acceleration of the 
form, Q(E p ) cx E~ 2 cxp(-E p / E max ) [s" 1 GeV" 1 ]. In or- 
der to calculate the spectra of cosmic ray protons and neutri- 
nos escaping from AGN jets we first need to obtain the ambi- 
ent spectrum of protons, N p (E p ) [GeV -1 ], in the shocked re- 
gions where the soft photon target density is sufficiently high 



for pion photoproduction to take place. From the ambient pro- 
ton spectrum we can then obtain the spectra of neutrinos and 
neutrons produced, Q n (E n ) and Q V (E V ) [s _1 GeV -1 ], re- 
spectively. Then taking account of the optical depth of the 
emission region for neutrons escaping photo-induced n — > p 
conversions, one can obtain from the neutron production spec- 
trum, the spectrum of cosmic ray protons resulting from the 
/3-decay of the neutrons escaping from the jet, Q a (E p ). 

The spectrum of protons on acceleration and the ambient 
proton spectrum are related by the proton loss time scale, 
tp{Ep), 

N p (E p ) ~ Q p {E p )t p (E p ) . (5) 

The processes mainly contributing to losses of protons are in- 
teractions with radiation, advection away from the shock re- 
gion of dimension R, and adiabatic energy losses if the emis- 
sion region expands. 

The advection time scale is expected to be i ac i v ~ R/((3 s hc) 
where /3 s h is the shock velocity and R is the dimension of the 
jet in the shocked region. In relativistic flows streaming away 
from a central source, this energy independent time scale is 
usually in competition with adiabatic energy losses of the pro- 
tons due to the expansion of the flow [16]. In a relativistic out- 
flow, characterized by a bulk Lorentz factor T and an open- 
ing angle 9, the expansion velocity in the co-moving frame 
(in units of c) is /3 CX rts TO for 6 < T" 1 , and /3 CX rts 1 oth- 
erwise. This leads to adiabatic cooling on a time scale i ac i w 
i?/(/3 ex c). For example, in GRB one can generally assume 
that /3 ex — 1, and also observations of superluminal motions 
in AGN jets are consistent with 9 ~ T -1 [17], thus /3 CX ~ 
1. In the following, we shall assume that adiabatic losses are 
relevant with 0.3 < /3 0X < 1. This has the important con- 
sequence that the lifetime of protons in the jet (or outflow) is 
limited to about one crossing time. The time scale for diffu- 
sive escape of protons is usually much longer (except, maybe, 
near the maximum proton energy), thus protons with E p -C 
£ max can be assumed to be confined in the emitter and do not 
contribute to the cosmic ray emission. 

The photon target spectrum will be assumed to have a 
power law shape, n(e) oc e - " -1 , extending to energies suffi- 
ciently above the threshold for photopion production by pro- 
tons of energy E p . Then, for a > the time scale for en- 
ergy loss by photohadronic interactions is asymptotically of 
the form 

t P1 (E p ) cx E~ a , (6) 

where t pl is understood as including Bethe-Heitler and pion 
production losses, the cooling time for pion production will 
be called i Pi7r > t P1 . For very flat target spectra, as for exam- 
ple in GRB at ultra-high proton energies, the photoproduction 
time scale is approximately constant [16], thus, Eq. (6) ap- 
plies with a = 0. We shall confine the discussion to the val- 
ues a = 1, relevant in AGN, and a = hereafter. Hence, for 
a = 1 we obtain 

t pl (E p ) = {E x /E p ){R/c), (7) 
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with tp^(E p ) w 1.5tp~ ( (E p ). Ei is the energy corresponding 
to unit optical depth for photoproduction losses, t p1 (Ei) = 
R/[ctp 1 (E 1 )} = 1. Setting i" 1 = t"^ + t£ we obtain, 

^V p (^p) = Qp(Ep)(R/c)[(E p /E 1 ) + a}- 1 (8) 
with a = max(/3 cx , /3 s h), and for Q p (E p ) oc i?~ 2 we have 



However, because the two break energies, aE\ and bE\, are 
probably very close, a ~ 6, we shall adopt a single break 
energy £^ = &Ei, and use the following approximations 



Qn(E ni L p ) oc L p exp 



-En 



E„ 



E n lE > 

E- 2 



(E n < E b ) 
(E b < E n ) 



(13) 



E- 2 (E p < aEi) 



n p (Ep)oc^ EiE j> 3 ^ >aE - 



(9) 



(Note that for clarity, here and in the next section we omit the 
exponential cut-off in the spectrum at E max .) Obviously, for 
a = the optical depth for photoproduction is constant, and 
Np(Ep) oc Qp(Ep). One can show that for typical photon 
densities in GRB fireballs r p7 < 1, and that Bethe-Heitler 
losses are unimportant, viz., t Pi7I w t pl [18,16]. 



Qcr(E p ,Lp) oc L p exp 



-E„ 



E n 



E7 1 



Ep 3 E b 



(E p < E b 
(E b < E p 



(14) 



Note that we have now put in explicitly the cut-off in the ac- 
celerated proton spectrum, and the proportionality with the 
proton luminosity of the source, L p . In Sect. IV, we shall 
relate L p to the to observed photon luminosities for spe- 
cific models of neutrino emission by AGN jets. Obviously, 
C. Generic cosmic ray proton and neutrino production spectra in GRB we simply have Q CT (E n ,L p ) « Q n (E n ,L p ) oc 

LpE~ 2 e-xp(—Ep/E mayL ). 

The production spectrum of muon-neutrinos will have the 
same broken power-law form as the neutron production spec- 
trum, and is related to it by 



The time scale for photohadronic production of neutrons is 



t 



P7 — >n ~ tp -7T K p / (N n ) 

E^ 1 , while for a 



For a = 1 
we have t 



this is t 



P7— m 



P7 

2.5t 



in- 



. 5ip^y oc 
This imme- 



diately gives the production spectrum of neutrons, 



(En) 

{E n < aEi) 
{E n > aE^ 



oc 



I E lE - 2 



Q^(E) = l-^l_Q n (E/r 1m ) 



(15) 



(10) 



Neutrons may escape to become cosmic ray protons. How- 
ever, because neutrons themselves suffer pion photoproduc- 
tion losses, the cosmic ray production spectrum will differ 
from Q n (E p ) above the energy, bE\ at which the optical depth 
for neutron escape, T n7 , is one. Neutrons can be considered as 
"absorbed" after they are converted into a proton, or after they 
have lost most of their energy in nj interactions, whichever 
time scale is shorter. For a = 1 this means r„ 7 ~ 2r p7 , giv- 
ing b w 0.5 for AGN jets, while for a = and typical GRB 
photon densities, t„ 7 w t P7 < 1, which means that neutron 
absorption is unimportant in GRB. 

We note that in a homogeneous spherical medium of radius 
R the optical depth decreases radially oc (1 — r/R) from its 
central value t to zero at r = R giving rise to the geometrical 
escape probability of an interacting particle propagating in a 
straight line 



where we count and together, and the correspond- 
ing spectrum of electron neutrinos at the source would be 
Qv e (E) w ^Q Ufi (E). Putting in the numbers given in 
Sect. II A, we find 



Q Vti (E) w 83.3Q„(25£) for a = 1 
Q„„ (E) w 416.Q„(25£) for a = 



(16) 
(17) 



V csc (t) : 



i(l-e- 
( 1 



r 



-l 



')/r 

T < 1 
T > 1 ' 



(ID 



resulting in the cosmic ray proton production spectrum being 
steepened above bE\, compared with Q n (E p ). For a = 1, the 
cosmic ray proton production spectrum is therefore 



Q C r(E p ) OC 



'Ep'E- 



E P - 2 
bEiE~ 



{E p < o£7i) 

(a£i < E p < bE{) 

(bE x < E p ) 



(12) 



We shall refer to Eq. (14) and Eq. (15) as the generic cos- 
mic ray and neutrino production spectra. We emphasize the 
strong dependence of the number of produced neutrinos per 
produced neutron on the assumed target photon spectral in- 
dex: at ultra-high energies, GRB produce about 5 times more 
neutrinos per neutron than AGN. We shall return to the impli- 
cations of this result at the end of the paper. 



III. PROPAGATION OF NEUTRINOS, PHOTONS, AND 
PROTONS OVER COSMOLOGICAL DISTANCES 

In this section, we discuss propagation of cosmic rays and 
neutrinos in an expanding Universe filled with the cosmic 
microwave background radiation. To illustrate the problem, 
we compare of energy-loss horizons of protons and neutri- 
nos. We shall then briefly discuss the physical problems con- 
nected to several approaches to cosmic ray propagation calcu- 
lations. Using the numerical propagation code described by 
Protheroe and Johnson [19], we then calculate the observable 
neutrino and cosmic ray spectra from a cosmological distri- 
bution of generic photohadronic sources, as described in the 
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last section. Here we assume that the sources are transpar- 
ent to neutrons, while protons are confined, and that gamma- 
rays are reprocessed in synchrotron-pair cascades until emit- 
ted in the energy range of 10 McV— 30 GcV. Using an extrap- 
olated cosmic ray spectrum which is consistent with present 
observational limits on the light component of cosmic rays 
(i.e. protons) as an upper limit on the possible extragalactic 
proton contribution, and the diffuse extragalactic gamma ray 
background (EGRB) observed by EGRET as an upper limit 
on the hadronic extragalactic gamma ray flux, we determine 
an energy dependent upper bound on the neutrino flux from 
cosmic ray sources with the assumed properties. Our result is 
compared with the energy independent bound on extragalactic 
neutrino fluxes recently proposed by Waxman & Bahcall [ 1 ] . 



A. Comparison of energy-loss horizons 

We wish to compare the distances that neutrinos, photons 
with energies below threshold for cascading in background 
radiation fields, and protons will travel through the Universe 
without significant energy losses. We define the energy loss 
horizon by 



X = cE/\dE/dt\, 



(18) 



such that for linear processes the energy is reduced to 1/e of 
its initial value on traversing a distance A. 

For gamma-rays below ^30 GcV and neutrinos, the 
energy-loss process is due to expansion of the Universe [20]. 
For simplicity, we adopt an Einstein-de Sitter cosmology (i.e., 
A = and £1 = 1), so that the horizon A can be related to a 
redshift z by the redshift-distance relation 



(19) 



where H a = 50/iso kms' 1 Mpc -1 is the Hubble constant. 
Since we require the distance for which the energy is reduced 
by a factor e during propagation, we have (1 + z) = e giving 
the horizon for redshift losses, 



X z = 



2c 



(1 



-3/2 



)■ 



(20) 



This is also the horizon for neutrinos and gamma-rays below 
~30 GeV. Normalizing the neutrino horizon to the radius of 
the Einstein-de Sitter universe, 



3iJ c 
2c 



(21) 



we obtain A~ 



Ay — X z 



0.78. 



In addition to redshift losses, extragalactic cosmic rays suf- 
fer energy losses from photohadronic interactions with cos- 
mic backgrounds, mainly the microwave background, and 
this is the reason for the Greisen-Zatsepin-Kuzmin (GZK) 
cut-off expected for a cosmic ray spectrum originating from 
a cosmologically homogeneous source distribution [21,22]. 



Photo-pion production and Bethe-Heitler pair production gov- 
ern the energy loss in different energy regimes due to their 
very different threshold energies. The Bethe-Heitler pro- 
cess limits the propagation of protons with energies E p > 
2m p m e c 4 //cT mbr w 4xl0 9 GeV to X p ~ 1 Gpc, while 
pion production reduces the horizon for protons with E p > 
mptn^c 4 1 'kT mhl w 5xl0 n GeV to A p ~ lOMpc. Again, 
in units of the radius of the Einstein-de Sitter Universe, the 
energy-dependent horizon for protons can be written 

X P (E) = [l/X z + 1/X p .bh(E) + l/Xp^E)}- 1 , (22) 

where the components expressing redshift, Bethe-Heitler and 
pion production losses, can be written as 



X z ~ 0.78 



A P) bh(-E) 
Xp^(E) 



0.27/i 50 exp(0.31/.Eio) 
5xl0" 4 ^5ocxp(26.7/£;io) 



(23) 



with Eio = E p /10 10 GeV. The approximations for A Pi bh 
and X Pt7r fit the exact functions determined numerically in [19] 
and the exact interaction kinematics within ~ 10% up to E p <~ 
10 12 GeV. 

The different energy-loss horizons for gamma rays and neu- 
trinos, and protons strongly affect the relative intensities of 
their diffuse isotropic background fluxes. This is true in par- 
ticular for evolving source populations such as quasars, galax- 
ies, or GRBs (if they trace star formation activity), since here 
most of the energy is released at large redshifts. Cosmic rays 
above the ankle (E CI ~ 3xl0 9 GcV) originate only from 
sources with redshifts z < 0.27, while neutrinos and gamma 
rays originate from sources within z v = z 7 ~ 1.7. This will 
give rise to the neutrino intensity being enhanced relative to 
the protons because of their larger horizon, and because of the 
evolution of the sources (e.g. quasars) with cosmic time (red- 
shift). 

We may illustrate the problem as follows. The basic method 
of calculating the approximate present-day diffuse fluxes of 
neutrinos, gamma rays below ~30 GeV, and cosmic rays of 
photoproduction origin, would be to integrate the contribu- 
tions from sources at redshifts up to those corresponding to 
the respective energy-loss horizons. Assuming a constant 
source number per co-moving volume element for simplic- 
ity, the resulting fluxes are proportional to V c (X)/d^(X) <~ A, 
where U C (A) and di(A) are the cosmological co-moving vol- 
ume and the luminosity distance, respectively, corresponding 
to the horizon A. Assuming photohadronic production of neu- 
trinos and cosmic rays in, e.g., a cosmological (non-evolving) 
distribution of AGN, the relative flux of neutrons (assum- 
ing no absorption) with energy E n , and corresponding neu- 
trinos with energy E v = 0.04£"„, is at the source given by 
[ElN v (E v )]/[E 2 n N n {E n )} = » 0.2. The same ra- 

tio will be observed in the integrated fluxes, as long E n < 
4xl0 9 GeV. For higher cosmic ray energies, the flux ratio 
must be multiplied by a factor A„/A p , which yields a flux ratio 
of -0.6 for E n ~ 10 10 GeV, and ~30 for E n ~ 10 11 GeV. 
Obviously, the differences would be much larger if we had 
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assumed strong source evolution which enhances the contri- 
bution from large distances. If we want to determine a neu- 
trino spectrum from an observed, correlated cosmic ray spec- 
trum from the same sources, the result must therefore approx- 
imately reflect the changes in the ratio of the energy loss hori- 
zons. 



B. Exact calculation of present-day neutrino and cosmic ray 
spectra 

The example above, of simply integrating the cosmic 
ray and neutrino contribution of cosmologically distributed 
sources up to the energy-loss horizon, disregards several im- 
portant aspects of particle propagation. Firstly, the particle 
number is not conserved in this method. Secondly, the en- 
ergy evolution of the particles is neglected, which removes 
the dependence on spectral properties of the source. Both 
caveats are removed in a method known as the continuous- 
loss approximation, which follows the particle energy along 
fixed trajectories as a function of cosmological distance (or 
redshift) [23,6]. For example, the trajectories for neutrino en- 
ergies would be simply E(z) — E (l + z). This method 
is exact for adiabatic losses due to the expansion of the Uni- 
verse (i.e., particle redshift), and is still a very good approx- 
imation for Bethe-Heitler losses, but it gives only poor re- 
sults for photo-pion losses. The reason for the latter is the 
large mean free path, and the large inelasticity of this pro- 
cess, which results in strong fluctuations of the particle en- 
ergy around its mean trajectory [24,25]. Cosmic ray transport 
in the regime where photo-pion losses are relevant is therefore 
best described by numerical approaches, either solving the 
exact transport equation [24], or by Monte-Carlo simulations 
[26]. Another important aspect concerning relative fluxes of 
cosmic rays, gamma rays and neutrinos is the fact that the in- 
teraction of the cosmic rays with the cosmic background ra- 
diation themselves produces secondary particles, which have 
to be considered as an additional contribution to the primary 
neutrinos and gamma rays. The full problem is treated by the 
cascade propagation code which has been described in detail 
by Protheroe and Johnson [19], and which we shall use also 
here. The code is based on the matrix-doubling technique for 
cascade propagation developed by Protheroe and Stanev [27]. 

For an input spectrum (dP ga i/dV c ){Q(E, z)) per unit co- 
moving volume per unit energy per unit time, the intensity at 
Earth at energy E is given by 

1(E) oc (24) 

Zmin 

where cLl and V c are luminosity distance and co-moving vol- 
ume, and M(E, z) are "modification factors" for injection of 
protons at redshift z as defined by Rachen and Biermann [6]; 
for neutrinos, M(E, z) = 1. The modification factors for pro- 
tons depend on the input spectra, and are calculated numeri- 
cally using the matrix method [19]. 
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FIG. 1. The observed all-particle cosmic ray spectrum spec- 
trum taken from the article by T.K. Gaisser and T. Stanev in the 
1998 Review of Particle Properties [81], and supplemented by Fly's 
Eye monocular data [28] (open circles at high energy), and recent 
AG AS A data [30] (rilled circles at high energy). Also shown are esti- 
mates of the cosmic ray proton component: based on the proton frac- 
tion estimates by [34] (thick lines with thick error bars, extended by 
thin lines which indicate the systematic error due to normalization to 
the all-particle data); Norikura data [82] (filled circles with large er- 
ror bars at 3 x 10 6 — 10 7 GeV); proton fraction estimated from KAS- 
CADE data [33] normalized to all-particle data (hatched band from 
10 6 — 10 7 GeV). The spectrum we adopt for the proton component 
which forms the upper bound to any extragalactic cosmic ray proton 
spectrum [cf. Eq. (25)] is shown by the dotted line. 

C. An abstract bound on extragalactic neutrino fluxes from 
neutron-transparent sources 

From the above considerations, it is obvious that one can 
use the observed cosmic ray spectrum to construct a correlated 
neutrino spectrum, under the assumption that all observed ex- 
tragalactic cosmic rays are due to neutrons ejected from the 
same sources as the neutrinos. A difficulty here is that the 
contribution of extragalactic cosmic rays to the total observed 
cosmic ray flux is unknown. Since neutrons convert into cos- 
mic ray protons, we can clearly consider only the proton com- 
ponent at all energies. Above the ankle in the cosmic ray spec- 
trum at E a ~ 3 x 10 9 GeV, observations are generally consis- 
tent with a "light" chemical composition, i.e., a 100% proton 
composition is possible. Since protons at these energies can- 
not be confined in the magnetic field of the Galaxy, they are 
also likely to be of extragalactic origin. At extremely high en- 
ergies, however, there is the problem that the event statistics 
is very low, and different experiments disagree on the mean 
cosmic ray flux at ^lO 11 GeV by one order of magnitude (see 
Bird et al. [28] for a comparison of the results until 1994 of 
the four major experiments Akeno, Havarah Park, Fly's Eye, 
and Yakutsk). This energy region is very important, since we 
expect here the existence of the GZK cutoff due to photopro- 
duction losses in the microwave background. Currently, no 
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FIG. 2. Spectra of (a) cosmic rays and (b) neutrinos after propagation through the Universe of input spectra for optically thin pion pho- 



toproduction sources for E u 



10 6 , 3 x 10 7 , 10 9 , 3 x 10 10 , 10 12 , and 3x 10 13 GeV, assuming galaxy evolution as described in the text. 



Spectra are normalized such that the cosmic ray intensity does not exceed the cosmic ray proton spectrum estimated from observations (dot- 
ted line in part a) and such that the neutrino energy flux does not exceed 0.5 of the observed photon energy density above 3 MeV. The dotted 
curve in part b joins the peaks in the neutrino spectra and forms our neutrino upper bound for optically thin pion photoproduction sources. The 
dashed line is the bound obtained by Waxman and Bahcall [1]. 



clear evidence for the existence of this cutoff has been found, 
and the results of at least two major experiments (Havarah 
Park [29], and AGASA 1998 [30]) are consistent with the re- 
sult obtained from a superposition of all experiments using 
a maximum likelihood technique [31,32], that is, a continua- 
tion of the cosmic ray spectrum as a power law ocE 2 ' 75 up to 
>3xlO n GeV (see also Fig. 1). 

The situation is even more difficult at lower energies: cos- 
mic rays are here assumed to be mainly of Galactic origin, and 
there is evidence that a considerable, maybe dominant frac- 
tion consists of heavy nuclei rather than protons. Around the 
knee or the cosmic ray spectrum at E cr ~ 10 6 -10 7 GeV, re- 
cent results from the KASCADE air shower experiment sug- 
gest that the fraction of heavy nuclei in the cosmic ray flux is 
at least ^30%, and further increasing with energy [33]. Also 
below the ankle, in the energy range 10 8 — 10 9 GeV, the anal- 
ysis of air shower data has produced tentative evidence of a 
composition change from heavy to light (with increasing en- 
ergy), supporting a dominantly heavy composition of cosmic 
rays between the knee and the ankle [34]. (Note that this re- 
sult is under dispute, and it has been shown that it depends on 
the Monte-Carlo simulation codes used to construct the air- 
shower properties in dependence of the primary particle mass 
[35]. These simulation codes involve particle interaction mod- 
els based on extrapolations many orders of magnitude above 
the energy range currently accessible with particle accelera- 
tors [36].) 

Using all the available data, we find that an extragalactic 
cosmic ray spectrum of the form 

N P}0bs {E) = 0.8 x (£7/1 GeV)- 2 - 75 cm- 2 s" 1 sr" 1 GeV" 1 
(3 x 10 6 GeV < E < 10 12 GeV) (25) 



is consistent with all data and limits on the cosmic ray proton 
flux (Fig. 1). It represents the current experimental upper limit 
on the extragalactic cosmic ray proton flux, which we shall 
use to construct an upper limit on the possible, diffuse extra- 
galactic neutrino flux. If it can be shown that the intensity of 
protons at 10 9 GeV and lower energies is below that assumed 
in this paper (dotted line in Fig. 1), then the neutrino bound we 
have constructed below 10 7 GeV would need to be reduced. 

To construct the neutrino bound, we assume test spec- 
tra of the form Q CI (E) = Q n (E) oc E- 1 cxp(- E / E max ) 
with 10 6 GeV<-E max <10 12 GeV. The corresponding neu- 
trino spectra are determined using Eq. (16). We assume 
a source distribution following the cosmological evolution 
function found for galaxies and AGN ( [37], see next sec- 
tion). For a given E max , the total contribution of cosmic rays 
and neutrinos is calculated using Eq. (24). The resulting spec- 
trum is then normalized so that its maximum reaches the cos- 
mic ray flux given by Eq. (25). By varying E max between 
10 6 GeV and 10 12 GeV, we then obtain the desired maximum 
flux of neutrinos consistent with the present cosmic ray data 
(see Fig. 2). We also consider the correlated gamma ray out- 
put, assumed to be twice the neutrino energy flux, and check 
it does not exceed the observed power-law component of the 
diffuse gamma-ray background (we estimate the background 
between 3 MeV and 30 GeV to be wl.5xl0" 5 GeV cm~ 2 
s _1 sr _1 , [38]). Note that the cosmic ray curve for E max — 
10 6 GeV does not reach the estimated cosmic ray proton spec- 
trum in Fig. 2(a) in order to avoid over-producing diffuse 
gamma rays. Note also that the propagated cosmic ray spec- 
tra cut off at or below 10 11 GeV for all E max values, and so 
our bound is insensitive to the assumed extragalactic cosmic 
ray spectrum above 10 11 GeV. 

The upper bound on the neutrino flux according to this 
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FIG. 3. Muon neutrino upper bounds for optically thin pion pho- 
toproduction sources (curve labeled t„ 7 < 1) and optically thick 
pion photoproduction sources (curve labeled t„ 7 > 1); the hatched 
range between the two curves can be considered the allowed region 
for upper bounds for sources with t„ 7 > 1. For comparison we 
show the bound obtained by Waxman and Bahcall [1] (for an evolv- 
ing source distribution). Predictions for optically thin photoproduc- 
tion sources are also shown: proton-blazar (Mannheim 1995 [10], 
Model A) - dotted curve, and GRB sources (Waxman and Bahcall 
1997, [18]) - dashed curve. Also shown is an observational upper 
limit from Frejus [83] and the atmospheric background [84]. 

model is given by the minimum of the cosmic ray and the 
gamma ray bound, shown in Figure 3 by the curve "t„ 7 < 1". 
We see that the cosmic ray limit starts to dominate the bound 
for E v > 100 TeV, which then decreases to a minimum at 
E v w 10 9 GcV, after which it rises again. This rise is a conse- 
quence of the strongly increasing ratio of neutrino and proton 
horizons above this energy, in conjunction with the assump- 
tions that the distribution of sources is homogeneous in space. 
A more realistic scenario might be that the lack of evidence 
for the GZK cutoff in the cosmic ray data is due to the domi- 
nant contribution of a local cosmic ray source. In this case, our 
estimate that [E 2 N v (E v )]/[E 2 N n (E n )\ oc A„/A p would not 
apply. Obviously, if an extended data set would confirm the 
existence of the GZK cutoff in the UHECR spectrum, the neu- 
trino bound would remain at the level of the Waxman-Bahcall 
result (also shown in Fig. 3, see discussion below) for E v > 
10 9 GeV. 

By assumption, the neutrino bound constructed this way ap- 
plies only to sources which are transparent to neutrons. For 
the opposite extreme, i.e., sources with a very high neutron 
opacity, t„ 7 » 1, it is still possible to set an upper limit us- 
ing the observed EGRB, assuming that the dominant part of 
the emitted gamma-radiation is in the EGRET range. This is 
shown by the line labeled "t„ 7 3> 1" in Figure 3. The range 
in between can be regarded the "allowed range" for the neu- 
trino emission from sources with r„ 7 > 1. In the next sec- 
tion, we will estimate specific neutrino spectra for AGN mod- 
els which imply that T ni {E n ) > 1 for high neutron energies. 



Our result may be compared with the cosmic ray bound to 
extragalactic neutrino fluxes recently proposed by Waxman & 
Bahcall [1]. Their bound was constructed using a different ap- 
proach: Waxman & Bahcall assume an E~ 2 input spectrum of 
extragalactic cosmic rays, and normalize the propagated spec- 
trum to the observed flux at E CI = 10 10 GeV. Consequently, 
their bound (determined for a source evolution oc (1 + z) 3 5 ) 
agrees with the one derived in this work at E v ~ 5 x 10 8 GeV, 
where the cosmic ray limit is most restrictive. 

We have chosen E^ 1 trial spectra with variable exponential 
cutoffs in order to be able to mimic the effect of the superpo- 
sition of spectra from various source classes. The pronounced 
peak in the energy flux associated with the trial spectra allows 
one to normalize the neutrino flux consistent with the exper- 
imental upper limit on extragalactic protons (oc E~ 2 - 75 ) at 
any chosen energy. Other hard spectra or delta-function dis- 
tributions as trial spectra would have yielded practically the 
same result. Although canonical AGN jets, which we dis- 
cuss in the next section, are an example for sources with E~ x 
(or similar) spectra, our result does not imply that we as- 
sume AGN jets to actually saturate the upper limit in gen- 
eral. As a matter of fact, at neutrino energies >10 9 GeV, a 
class of photohadronic sources saturating our bound has nei- 
ther been suggested nor does its existence seem likely on the 
basis of current knowledge. By restricting their source spec- 
tra to E~ 2 , Waxman & Bahcall have constructed a bound for 
neutron-transparent sources which is probably closer to cur- 
rent models for cosmic ray and neutrino production than our 
general upper limit above 10 19 eV. Examples for such mod- 
els are the model for diffuse neutrino fluxes from GRBs pro- 
posed by the same authors [18,1] (note that this prediction as- 
sumes no cosmological evolution for GRBs), and Model A 
in Mannheim [10] which was constructed using the cosmic 
ray limit with the assumption that the emerging neutrons at 
~10 10 GeV contribute to the extragalactic cosmic ray spec- 
trum (both shown in Fig. 3). 

IV. DIFFUSE NEUTRINO SPECTRA FROM AGN JETS 

In this section, we shall consider the spectra of cosmic ray 
protons and neutrinos emerging from jets of two classes of 
gamma-ray emitting AGN, i.e., BL Lac objects and radio- 
quasars, which are usually combined as the class of blazars. 
We shall use the estimated optical depths of gamma rays to 
photon-photon pair production in a typical AGN of each type, 
to infer the corresponding neutron-photon optical depths in 
these objects. We shall show that high luminosity AGN (like 
3C 279) can be expected to be opaque to neutrons at energies 
above about 10 8 — 10 9 GcV, while low luminosity BL Lacs 
(like Mrk501) must be transparent to neutrons at all ener- 
gies. Then assuming a model for the luminosity dependence 
of the optical depths, and the local luminosity functions of 
BL Lacs and quasars, we shall estimate the form of the pro- 
duction spectra of cosmic ray protons and neutrinos per unit 
volume of the local Universe. Applying the source evolution 
functions found for BL Lacs and radio-quasars, respectively, 
we shall derive model estimates for the diffuse neutrino con- 
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tribution from these sources which are compatible with cos- 
mic ray limits. We shall also construct an upper bound for 
the contribution of AGN jets, using the same method as in the 
previous section, but for the appropriate generic source spec- 
tra, Eqs. (14) and (15), varying the break energy Eb over the 
range allowed by the models. 

A. Cosmic ray proton and neutrino production spectra from 
blazars 

As our starting point, we assume a target photon spec- 
trum with index a = 1 which we have already seen leads 
to T P j(E p ) oc E p . A similar energy dependence applies to 
gamma rays interacting with the same photons by photon- 
photon pair production (77 — > e + e~), and so for a = 1 we 
have 



ipj(Ep) 



2m e c 2 E p 



(ct 77 } 77 V win [m p + m w /2] 
5xl(r 4 T 77 ([4xlCr 6 ]£g . 



(26) 



Here we have used (K p a P7 ) / (c 77 ) s=s 300/ibarn/(7T from 
averaging over a photon spectrum with a = 1, where (Kp(T P7 ) 
includes Bethe-Heitler pair production, cr 77 is the total cross 
section for the process 77 — > e+e~ [39], and or is the Thom- 
son cross section. Using T ni (E n ) w 2r ln (E n ), and assuming 
that T ni (E) oc E holds for a range 10 5 £" 7 < E n < £ ma x, we 
obtain the relation 



,(£„)/t 77 (£ 7 ) m 4 x 10- 9 E„/E 7 



(27) 



We apply this relation to two reference AGN: the BL Lac ob- 
ject Mrk501, and the quasar 3C279. The combination of 
the observed TeV spectrum and the EGRET flux limits for 
Mrk501 [40,41] gives rise to our assumption of a break en- 
ergy at about . 3 — 1 TeV in the context of photohadronic mod- 
els, while for 3C 279 the EGRET spectrum allows for a break 
at ^3— 10 GcV [42]. Note that internal opacity due to the 
presence of low-energy synchrotron photons with spectral in- 
dex a, which plays a crucial role in photohadronic models for 
the gamma ray emission, gives rise to a spectral steepening by 
E~ a , and not to an exponential cutoff, above the energy where 
t 77 = 1, see Eq. (11). A TeV power law being steeper than 
the gamma ray spectrum at EGRET energies could thus imply 
optical thickness in spite of reaching some 25 TeV as seen by 
the HEGRA telescopes [43]. Clearly, lower values for the op- 
tical depth are obtained if the data are interpreted with models 
which assume that the observed steepening is due to the radi- 
ation process itself (e.g. synchrotron-self-Compton emission) 
[40,41]. Applying Eq. (27) we find that the break energy in 
the cosmic ray production spectra should be Eb ~ 10 11 GeV 
in Mrk501, and E b ~ 10 9 GcV in 3C279. Thus, UHE cos- 
mic rays can escape from Mrk501 under optically thin condi- 
tions, and our upper bound is not affected by the actual value 
of the gamma ray optical depth in this source. 

We note that the break energies derived above depend on 
the assumption of an undistorted power law target spectrum of 
photons. This assumption might be reasonable for BL Lacs, 



which show power law photon spectra extending from the in- 
frared (relevant for pj interactions) into the X-ray regime (rel- 
evant for 77 absorption of ~10— 100 GeV photons, assuming 
a Doppler boosting of the emission by a factor S ~ 10). In 
high luminosity radio-quasars, however, this is not the case 
(see, e.g., [44]), and relation (27) does not necessarily hold. 
Moreover, in these sources the dominant target spectrum for 
pj and 77 interactions could be given by the external ac- 
cretion disk photons, forming roughly an (unboosted) power 
law spectrum with a = 1 up to ~10 eV [45], where it drops 
by about one order of magnitude (e.g. [46]). This still im- 
plies a relation like Eq. (27), but with r„ 7 increased relative 
to t 77 by a factor of 10 — consequently, the possible break 
in the gamma-ray spectrum of 3C 279 would correspond to 
Eb ~ 10 8 GeV. A principal lower limit to the break energy 
in EGRET sources is set at Eb ~ 10 7 GeV, since otherwise 
EGRET photons (>1 GcV) could not be emitted. 

Using cosmic ray proton production spectra with a break 
at ^10 8 GcV (see Fig. 4) will, of course, have a strong ef- 
fect on the neutrino bound implied by cosmic ray data. Since 
the cosmic ray proton spectrum of the source above the break 
drops faster than the upper limit spectrum on cosmic ray pro- 
tons, Eq. (25), the bound at a neutrino energy E v > Eb /2b is 
essentially set by the cosmic ray flux at the break energy Eb, 
rather than by the more restrictive flux at 2bE v . If all sources 
had the same Eb, the bound would be increased roughly by 
a factor (2bE u / Eb) ' 75 . Clearly, the assumption of optical 
depths allowing break energies of 10 8 GeV or below in lu- 
minous quasars is not directly supported, but only consistent 
with current observations. The maximal neutrino intensity for 
AGN jets, which we shall derive below on the basis of this as- 
sumption, is therefore at the currently allowed maximum for 
the adopted source model, and likely to be lowered when more 
detailed gamma-ray data become available. 

As a general limitation of this approach it must be noted 
that most AGN jets have not been detected by EGRET. If 
their lack of high energy emission is due to a large gamma- 
ray opacity, this would also imply very low break energies for 
their ejected cosmic ray spectrum. One such possible class 
of sources are the gigahertz -peaked sources (GPS) and com- 
pact steep-spectrum (CSS) quasars [47], which make up about 
40% of the bright radio source population. Objects such as 
these could thus produce a diffuse neutrino intensity at the 
level of the EGRB or maybe even above, without violating 
any constraint. Also low-opacity sources for which the max- 
imum proton energy remains much less than its theoretically 
allowed value and, in particular, below the value of the break 
energy, could lead to a higher neutrino intensity at low E v 
(limited only by the EGRB). We did not consider this possi- 
bility by keeping E Ptinax fixed. 



B. Blazar luminosity functions and generic models for the 
neutrino contribution from BL Lacs and radio-quasars 

In order to obtain a parametrization of blazar neutrino spec- 
tra, we need to express T n7 as a function of the blazar lumi- 
nosity L, given at some frequency. Here we take into consid- 
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FIG. 4. Form of the spectra of neutrinos and cosmic rays escap- 
ing from optically thick photoproduction sources, Eqs. (14) and (15). 
The example shown here has a break energy of 10 8 GeV, and an ex- 
ponential cut-off above 10 11 GeV. 



eration that blazars are assumed to be beamed emitters, with a 
Doppler factor S <~ 10. The optical depth intrinsic to the emis- 
sion region, is proportional to L t / R, where L t is the intrinsic 
target photon luminosity and R is the size of the emitter. Since 
blazars are strongly variable objects [48], we can use the vari- 
ability time scale, T var , to estimate the intrinsic size by R ~ 
T val cS. Using also the relation between intrinsic and observed 

Al- 



luminosity, L = L t S , we obtain r, 



717 



LT-JS~ 5 



though there is no detailed study of a possible systematic de- 
pendence of T var on L, the observations are compatible with 
no such correlation existing — for example, variability time 
scales of order 1 day are common in both, moderately bright 
BL Lacs like Mrk 501 or Mrk 421, and powerful quasars like 
3C279, PKS 0528+134, or PKS 1622-297, whose optical lu- 
minosities differ by at least three orders of magnitude [48]. 
For the Doppler factors, there is no evidence for a systematic 
dependence on L either [17], although unification models for 
blazars and radio galaxies (see next section) suggest that BL 
Lacs are on average slightly less beamed ({5} ~ 7) than radio- 
quasars ((5) ~ 11) [49]. Therefore, we may assume that for 
both BL Lacs and quasars, r„ 7 (_E, L) oc L holds on average, 
and that T nj (E, L, BL Lac) <~ 10t„ 7 (_E, L, quasar). It is in- 
teresting to note that this relation would imply a relation of the 
break energies of Mrk501 (L opt ~ 10 44 erg s" 1 ) and 3C 279 
(L opt ~ 10 47 erg/s) as£ 6 (3C279) - 10~ 2 £; fc (Mrk 501), 
consistent with our estimate obtained for intrinsic absorption 
from the gamma ray spectral break. Of course, this does not 
rule out the possibility that Eb might be systematically lower 
in 3C 279 due to external photons, as argued above. 

To determine the contribution of all blazars in the Universe, 
we have to relate the proton luminosity L p to the blazar lu- 
minosity L in some frequency range, and then integrate over 
the luminosity function, dN / dL, determined for the same fre- 
quency range. Here we have to distinguish between the lumi- 
nosities in the energy range where the target photons are, L , 



and the luminosity of the gamma-rays, L 7 , which are here as- 
sumed to be produced by hadronic interactions. Obviously, 
the latter implies L p oc L 7 , but on the other hand r„ 7 oc 
L Q . We also have to distinguish between the two classes of 
blazars: while for BL Lacs L 1 oc L , observations rather sug- 
gest that for quasars L 7 oc L 2 D [44]. This leads to the follow- 
ing simple models: 

For BL Lacs, we use the X-ray luminosity functions given 
by Wolter et al. [50] for X-ray selected BL Lacs, 



1.6 



gWbl / dLx oc L x 

(SxlO^ergs" 1 < L x < SxlO^ergs" 1 ), 



(28) 



and use the relations L p oc Lx 



and Eb oc L x with Eb = 
SxlO^ergs" 1 (Mrk501). Here we 



10 11 GeV for L x 
have assumed that L Q oc Lx- 

For quasars, we use the EGRET luminosity function given 
by Chiang & Mukherjee [51], 



-2.2 



dAfq/dLj ^ 7 

(10 46 ergs < L 7 < 10 48 



(29) 



ergs 



-1/2 



and the relations L p oc L 1 , and Eb oc L 7 1/z where we 
consider two possible normalizations for L 7 = 10 48 ergs _1 
(3C279), which are Eb = 10 9 GcV in case that jet-intrinsic 
photons dominate the target field, and Eb = 10 8 GeV for the 
assumption that external photons are the dominant target. For 
illustration, we show in Fig. 4 cosmic ray and neutrino spec- 
tra on emission for E b = 10 8 GeV and £ max = 10 11 GeV. 

Then we obtain the form of the production spectra of cos- 
mic ray protons and neutrinos averaged over the local universe 
due to these two classes of AGN, 



(QcrAE)) 



J Q cr ,u{E, L)(dJ\f/dL)dL 
jL(dAT/dL)dL 



(30) 



where the input spectra Q cr , Q v are given by Eqs. (14) and 
(15), with E max = 10 11 GeV, and Eb, L p given as functions 
of L as discussed above. 

To integrate properly over redshift, we note that while 
quasars show strong evolution similar to galaxies (see below), 
BL Lacs show little or no evolution [52], and we shall take 
this into account when propagating these spectra through the 
Universe from large redshifts. For quasars, their luminosity 
per co-moving volume has a pronounced peak at redshifts of 
z ~ 2, and declines or levels off at higher redshifts [53]. We 
shall assume that this effect is due to evolution of the num- 
ber of quasars with z, rather than evolution of the luminos- 
ity of individual sources, which keeps the production spectra 
(Qa,u(E)) independent of z. A particular parametrization of 
the redshift-dependence of the (co-moving frame) UV lumi- 
nosity density of AGNs as inferred by Boyle and Terlevich 
[37], assuming an Einstein-de Sitter cosmology and h 50 = 1, 
is given by 



dP^ 



dV c 



[(l + z)/2.9] 3 ' 4 (z<1.9) 
P { 1.0 (1.9 < z < 3) 

exp[-(z - 3)/1.099] (z > 3) 



(31) 
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where P = (3.0 ± 0.3)xl0 44 crg s" 1 Mpc~ 3 . Clearly, the 
normalization plays no role here since L p is adapted to match 
the cosmic ray flux at earth. So for BL Lacs (no evolution), 
we simply use dP ga i/dV c = 1. 

The result is shown in Figure 5. As expected, the diffuse 
neutrino fluxes from AGN jet models exceed the bound for op- 
tically thin sources, but fall into the allowed region for sources 
optically thick for neutron emission. The BL Lac contribu- 
tion falls below the bound, because it was derived for a non- 
evolving source distribution — we note that it still exceeds 
the corresponding Waxman-Bahcall bound for the case of no 
evolution. In addition to the models discussed above, we have 
also constructed and estimated an upper bound on the diffuse 
neutrino contribution on AGN jets. Here we used the same 
method as for the construction of Fig. 3, but implying input 
spectra of the form of Eqs. (14) and (15), with variable Eb 
and fixed -E max — 10 11 GeV. The break energy was then 
varied in the range 10 7 GeV < E b < E max = 10 11 GeV, 
and the normalization chosen such that the superposed spec- 
tra approximately represented the upper limit on the extra- 
galactic cosmic ray spectrum (Fig. 5, left panel). We note 
that this upper bound corresponds within a factor of 2 with 
the prediction of a previously published model by Protheroe 
[45]. Other published models, for example by Halzen & Zas 
[ 1 1 ] , or Mannheim [ 1 0, Model B] , exceed our bound by about 
one order of magnitude at energies E v > 10 8 GeV, but their 
predictions for the important energy range below ^10 7 GeV 
are compatible with this bound. We return to the discussion 
of these models later, after we discuss the effect of magnetic 
fields in the AGN environment in the next section. 



V. MAGNETIC FIELDS AND THEIR IMPACT ON COSMIC 
RAY PROPAGATION 

As discussed in Sect. II, protons accelerated inside AGN 
jets are likely to be trapped in the jet to be released later near 
the end of the jet, and will consequently suffer severe adia- 
batic losses as a result of jet expansion. The bound that we 
calculated for optically thin photoproduction sources shown in 
Fig. 3 may of course be exceeded for optically thick sources. 
However, it may even be exceeded for optically thin photo- 
production sources if the cosmic rays resulting from photo- 
produced neutrons are trapped, or suffer severe adiabatic de- 
celeration in the large-scale magnetized environment such as 
the host galaxy and its halo, a galaxy cluster, or a superclus- 
ter. We emphasize that although there is now useful informa- 
tion available about magnetic fields in galaxies and clusters, 
the magnetic field structure and topology is not sufficiently 
well known for us to predict reliably magnetic trapping and 
adiabatic losses in host galaxies and clusters. However, we 
shall discuss in this section the fate of the cosmic rays re- 
sulting from photoproduced neutrons, and show that in some 
plausible scenarios, our bound for optically thin photoproduc- 
tion sources can be exceeded for neutrinos from optically thin 
photoproduction sources. 

X-ray observations and measurements of extragalactic 
Faraday rotation suggest that structures surrounding compact 



AGN jets carry magnetic fields of the order of 0.1-10 /iG [54]. 
These can influence the propagation of cosmic ray protons in 
essentially two ways: (a) particles may be physically confined 
in the structure for a time t > in = l/-ff , or (b) the diffu- 
sive escape of the particles can lead to adiabatic energy losses. 
Obviously, magnetic fields on scales larger than the mean free 
path of a neutron for /3-decay, l n « 10 kpc(-E„/10 9 GeV), 
also affect cosmic rays which are ejected as neutrons from the 
source. Here we shall discuss what influence these effects can 
have on the strength of the measured cosmic ray flux relative 
to the corresponding neutrino flux. We shall estimate the crit- 
ical energy E* below which our bound could plausibly be ex- 
ceeded for optically thin photoproduction sources for a num- 
ber of scenarios, in particular for clusters of galaxies and ra- 
dio galaxies hosting AGN. 



A. Particle confinement in clusters and superclusters 

Clusters of galaxies have been recently discussed in the lit- 
erature as a possible "storage room" for cosmic rays [55-57], 
because of their relatively strong magnetic fields (Bq > 1 /iG) 
extending over large scales, i.e. cluster radii, R c \ > 1 Mpc. 
The time scale for diffusive escape of cosmic rays in a turbu- 
lent magnetic field of homogeneous strength B within a cen- 
tral radius R is given by 



, R o 

tesc ~ 



3eB Rl 



D cE cr \(E CT ) 



(32) 



where we have used for the diffusion coefficient D = ^Xr^c, 
and A (.E'er) is the scattering length in units of the Larmor ra- 
dius, tl = -E cr /e-Bo of a cosmic ray proton of energy E cr . 
The function A depends on the turbulence spectrum of the 
magnetic field, [8B(k)\ 2 oc k~ y , where k is the wavenumber. 
Cosmic ray scattering is dominated by field fluctuations on the 
scale of the Larmor radius, i.e., k ~ . This implies that 



a(e; c 



eB cio 



E r 



for -E cr < eBa a 



(33) 



where ao = fc~ in is the largest scale of the turbulence, i.e., the 
"cell-size" (or "reversal-scale") of the magnetic field. For k > 
clq 1 , the magnetic turbulence in clusters of galaxies seems to 
be well described by the Kolmogorov law for fully developed 
hydrodynamical turbulence, i.e., y = |. This can be easily 
seen from relating the typical turbulent magnetic field and cell 
size found from Faraday rotation measurements, ao ~ 20 kpc 
and Bo ~ 1 /iG [58,59], to the diffusion coefficient found 
for electrons of energy ~1 GeV from the synchrotron radio 
emission spectrum [60], D <~ 2xl0 29 cm 2 s~ x . For E„ > 
eBao, i.e. tl > ao, the particle motion is a random walk 
with scattering angles ^a / in each step. This can also be 
approximately described by Eqs. (32) and (33) by setting y = 
-1. 

Confinement of cosmic rays over the cluster radius, i? = 
i?ci ~ 3 Mpc, is then obtained for i CS c > tn, corresponding 
to a critical energy 
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FIG. 5. (a) Comparison of neutrino spectra for optically thick pion photoproduction sources with neutrino upper bounds obtained in the 
present work and by Waxman and Bahcall [1]. Neutrino intensities obtained in the present work plotted are: maximal superposition of spectra 
having Et in the range 10 7 — 10 11 GoV (thick solid curve); maximum source spectra averaged over the quasar luminosity function assuming 
t„ 7 oc L 1/2 with break energies corresponding to L = 10 48 erg s _1 at E b w 10 8 GeV (thick dashed curve) and at E b w 10 9 GeV (thick 
dotted curve); maximum source spectra averaged over the BL Lac luminosity function assuming r„ 7 oc L with the break energy corresponding 
to L = 3 x 10 44 erg s _1 at Eb ~ 10 11 GeV (thin dot-dashed curve). Also shown is the prediction by Protheroe [45] for an external photon 
optically thick proton blazar model normalized down in accordance with the recent estimates of the blazar contribution to the diffuse gamma 
ray background [53] (thick chain curve), and the bounds obtained by Waxman & Bahcall [1] with and without source evolution (note that the 
maximum BL Lac intensity was calculated assuming no evolution). Other symbols shown correspond to Figure 3. (b) Cosmic ray intensities for 
the neutrino intensities obtained in the present work shown in part a (key to these curves is as in part a). The upper bound to any extragalactic 
cosmic ray proton spectrum is shown by the thin line, which corresponds to the dotted line in Figure 1. 
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(34) 



provided ao > 1 kpc(i? c i/3 Mpc) 2 . This assumes that the 
magnetic field strength is homogeneous over the entire cluster. 
However, if it dropped to <~ O.lyuG, and ao ~ 200 kpc at the 
edge of the cluster then E^ iS would be a factor ^10 3 lower. 
On the other hand, some observations suggest also larger mag- 
netic fields on smaller reversal scales [61,62], which would 
imply higher confinement energies. 

The scenario above assumes that the background plasma 
filling the cluster is at rest. However, simulations of structure 
formation suggest that clusters of galaxies are accreting extra- 
galactic hot gas [63], forming inflows of typical speeds u acc > 
300 km s -1 downstream of an accretion shock near the outer 
radius of the cluster [7]. Since particles diffuse relative to the 
plasma flow, we require that in order to escape from the clus- 
ter, the average "speed" of particle diffusion, ^D/R c \, ex- 
ceeds the inflow velocity. If we assume the lower value for the 
magnetic field at R c \ <~ 3 Mpc, i.e. B ~ 0.1 ^G and a <~ 
200 kpc, we obtain a critical energy 



e: 



2xlO b GeVx 



(35) 



B(R ci y 




^cl ^acc 


3 


" a(R c i) ' 


0.1 /iG 




3 Mpc 300 kms" 1 . 




200 kpc 



Of course, with higher magnetic fields and a shorter field- 



reversal length a , E* cc would be higher. Similar processes 
would occur in the cooling flows observed in some rich clus- 
ters harboring powerful radio galaxies, but the confinement 
energies for typical cooling flow parameters [64,65] are even 
lower. Also the effect of drift in the nonhomogeneous cluster 
field [1] is small compared to diffusion, and can not lead to 
larger confinement energies. 

Ensslin et al. [57] point out that simple diffusion, as as- 
sumed above, may not be the best description of cosmic ray 
propagation in clusters of galaxies, and suggest that a one- 
dimensional random walk along static, but randomly tangled, 
magnetic field lines may be more realistic for particles with 
?*l <C a . In this case the field line topology may be con- 
sidered as arising from a three-dimensional random walk with 
steps of size a , yielding an effective diffusion length of R = 
R^/cio- For the escape time we then obtain 



2ca ! 



eB a 



E r 



1/3 



> t H for E cr < eB a , (36) 



which means for typical cluster parameters that essentially all 
cosmic rays with E cr < 3xl0 9 GeV are confined to the clus- 
ter. Of course, this result neglects cross-field diffusion so that 
the maximum confinement energy is probably lower, but it 
could still be higher than the result obtained in Eq. (34) for the 
case of simple diffusion, in particular in the likely case that 
the field strength decreases with the distance from the cluster 
core. 
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Confinement of cosmic rays in clusters would lead to a de- 
crease of the cosmic ray flux measured at earth relative to the 
corresponding neutrino flux, causing an increase of the neu- 
trino bound. It is important to note, however, that cosmic ray 
confinement may exist even on larger scales, i.e., superclus- 
ters. It has been shown that magnetic fields ~ 0.1 /iG in super- 
clusters are consistent with observations, and expected in sim- 
ulations of structure formation which also predict accretion of 
gas with speeds v acc ~ 1000 kms -1 [66,67]. The cosmic ray 
confinement energies for these larger scales (R ~ 10 Mpc, 
a ~ 1 Mpc) could therefore be >10 8 GeV, where again topo- 
logical aspects connected to the detailed structure of the field 
may allow even higher values. Since our Galaxy itself is lo- 
cated in a supercluster, this latter scenario would tend to in- 
crease the cosmic ray flux relative to the corresponding neu- 
trino flux, which fills the Universe homogeneously — thus it 
could actually decrease the bound below E v ~ 10 7 GeV. The 
net effect of confinement in clusters and our local superclus- 
ter is obviously strongly model dependent, and therefore dif- 
ficult to estimate. 

An effect of the unknown structure of magnetic fields on 
galaxy cluster and supercluster scales may therefore most 
likely affect the observed extragalactic cosmic ray flux at ener- 
gies below ~2xl0 6 GeV, corresponding to E v ~ 10 5 GeV. 
On the other hand, we can exclude an effect only for energies 
E CI > 3xl0 9 GeV {E v > 10 s GeV), but we note that cos- 
mic ray confinement at energies higher than ~ 10 s GeV re- 
quires extreme assumptions on the strength and topology of 
the magnetic fields. We therefore agree with the conclusion of 
Waxman & Bahcall [1] that at E v - 5xl0 8 GeV an effect of 
large scale magnetic fields on the relation of cosmic ray and 
neutrino fluxes cannot be expected. 



B. Adiabatic losses in expanding radio lobes and halos 

Although many galaxy clusters have powerful radio galax- 
ies at their centers, it is also a fact that most powerful radio 
galaxies are not found in such environments [68,69]. For a ra- 
dio galaxy located in the normal extragalactic medium, evi- 
dence has been found that a pressure equilibrium with the ex- 
ternal medium cannot be obtained within the lifetime of the 
source (<10 8 yr) [70]. Therefore, it can be concluded that 
the lobes must expand. For a sample of powerful double-lobe 
(or FR-II) radio galaxies, lobe propagation velocities of order 
10 4 kms -1 have been inferred [71]. Since the aspect ratio of 
the sources is found to be independent of their size, consis- 
tent with a propagation of the lobes along a constant opening 
angle of ^10°, expansion velocities of the lobes and the con- 
necting "bridges" of w ox ~ 1000 kms -1 can be inferred. 

Neutrons produced by pion photoproduction interactions at 
an acceleration site near the base of the jet will be beamed 
preferentially along the jet direction decaying farther out 
along the jet or in the radio lobe, where they decay. The result- 
ing cosmic ray protons will then be advected with the outflow- 
ing plasma on a time scale of the galaxy life-time ~10 8 yr <C 
in- Additionally, the protons perform a random walk in the 
magnetic field, which may even decrease their confinement 



time. Expanding radio lobes can therefore not confine cosmic 
rays indefinitely. However, in plasma outflows all particles 
which are isotropized in the flow due to scattering with plasma 
turbulence will experience adiabatic losses before their re- 
lease. Kinetic theory implies that the particle energy at ejec- 
tion from the flow is related to the injected particle energy by 



E\ n 



R n 



(37) 



Here, i? max is the radius of the outer termination shock of the 
flow, and i?i„j is in general given by some minimum radius 
i? m i n where the outflow starts. For the case of neutron ejection 



from a central source, as considered here, 



l n if l n ^> 



Rmin, and 



i? m i n otherwise. Note that the lobes we 



are considering here are much more extended than the jets in 
which the particles are accelerated. 

The observed synchrotron spectra from radio lobes imply 
typical magnetic fields in the range 10—50 /iG, turbulent with 
a maximum scale (or cell-size) of ^0.5 kpc [72,71]. The 
observed asymmetric depolarization in double radio galaxies 
leads to the suggestion that the magnetized plasma around the 
radio galaxy extends into a halo of radius ^300 kpc, with B ~ 
0.3 /xG and a cell size of ~5kpc at a radius R ~ 100 kpc [73]. 
The properties of the magnetic field in the central lobe and the 
halo can be connected by assuming that magnetic field and 
cell size scale as B = B (R/Ro)~ 2 and a = a (R/R ), re- 
spectively, with Bo ~ 30 /iG, ao ~ 0.5 kpc and Rq = 10 kpc. 

Within radius Ro, we assume the properties of the turbulent 
magnetic field are constant. This corresponds to the assump- 
tion of an isotropic magnetic field expanding in a plasma out- 
flow with i? m in = Ro and R max <~ 300 kpc, which will be 
used as a working hypothesis in the following. We can con- 
sider cosmic ray protons as nearly isotropized in the plasma if 
r L (i?) < a(R), corresponding to energies E cr (R) < i?* d (i?) 
with 



E* d (R) = 10 iU GeV x 



B 



30 mG 



a 



0.5 kpc 



R 



10 kpc 



(38) 



for 10 kpc < R < 300 kpc. To justify the assumption of 
near-isotropy below E* d (R), we have to show that advection 
in the flow indeed dominates over diffusion, i.e., that the cen- 
ters of the diffusive cloud of cosmic ray protons move ap- 
proximately steady with the flow. Indeed, for the parame- 
ters adopted above we find for the "speed" of diffusive escape 
D/R ~ ^a(R)v cc u/ R < v cx , where u co n < ^cis the average 
velocity of the particles to cross the cell, considering that also 
this motion is diffusive if the magnetic turbulence proceeds to 
smaller scales as expressed in Eq. (33). 

Consequently, we can consider the cosmic rays with E cr < 
i?* d (R) as are advected with the flow at a distance R from the 
center of the galaxy, so that they suffer adiabatic losses fol- 
lowing Eq. (37) leading to E cr (R) oc R -1 . This implies that 
E CI /E* d is independent of R, i.e. that cosmic rays confined in 
the outflow at some radius R remain confined for larger radii. 
If we consider cosmic rays which are ejected as neutrons, the 



13 



radius where they couple to the magnetic field is given by the 
/3-decay mean free path, l n « 10 kpc(£V,./10 9 GeV). The re- 
sulting protons are confined to the outflow and are subject to 
adiabatic losses provided E n < E* d {l n ), or 



E n < 3xlO y GcV 



r B 


1/2 


a 


1/2 


30 [iG 




0.5 kpc 





(39) 



The energy on ejection from the outflow is then E C j « E n /30 
for E n < 10 9 GeV, and E ci = ££/(3xl0 10 GcV) for 
10 9 GeV < E n < 3xl0 9 GeV. For E n > 3xl0 9 GeV, cos- 
mic rays (neutrons and protons) traverse the lobe/halo along 
almost straight paths and adiabatic losses do not apply. (Note 
that the corresponding energies where these modifications 
may affect the neutrino bound are E v ~ E n /20, thus E v < 
10 8 GeV.) 

Energy losses of cosmic rays in the lobes and halos of ra- 
dio galaxies are of particular relevance for models of neu- 
trino production in AGN jets, which we discussed in Sect. IV. 
These models apply to radio loud AGN, which are likely to 
be the beamed counterparts of radio galaxies [49]. The two 
classes of AGN discussed in the last section correspond to 
the two Fanaroff-Riley (FR, [74]) classes of radio galaxies: 
radio quasars might be associated to the powerful double- 
lobed FR-II radio galaxies, while BL Lac objects might cor- 
respond to the less luminous FR-I radio galaxies which gen- 
erally have diffuse lobes centered around the AGN. The pa- 
rameters used above for the lobes/halos of radio galaxies were 
mainly obtained from observations of FR-II radio galaxies or 
radio quasars. Therefore, the cosmic ray ejection from ra- 
dio quasar/FR-II sources can be expected to be diminished 
by more than an order of magnitude below ^10 9 GeV. For 
the less luminous FR-I galaxies, it could be that the magnetic 
fields, turbulence scales and halo sizes used above are overes- 
timated, so that E* d could be lower by about one or two or- 
ders of magnitude for these sources. 



VI. DISCUSSION AND CONCLUSIONS 

We have looked at the problem of to what extent a possi- 
ble flux of extragalactic neutrinos in the broad energy range 
10 5 GeV < E v < 10 11 GcV is bounded by the observed flux 
of cosmic rays and gamma rays. As the minimum contribu- 
tion to the cosmic rays we consider the neutrons produced 
together with other neutrals such as gamma rays and neutri- 
nos in photohadronic interactions. The most restrictive bound 
arises for sources which are transparent to the emission of 
neutrons and where the protons from the decaying neutrons 
are unaffected by large-scale magnetic fields. The bound for 
this case is approximately in agreement with the bound previ- 
ously computed by Waxman & Bahcall [1] in the neutrino en- 
ergy range 10 7 GeV ?SE V < 10 9 GeV but is higher at lower 
and at higher energies. The difference is a consequence of the 
different approaches: Waxman & Bahcall assume a particular 
model spectrum, which could explain the observed UHECR 
spectrum, but which has a pronounced GZK cutoff that is not 
clearly confirmed by the present low statistics data, and which 



does not reproduce the steep slope of the observed CR spec- 
trum below ^10 10 GeV. It thus leaves room for additional 
contributions from extragalactic sources (with different spec- 
tra) outside the narrow energy range where it matches with the 
observed cosmic ray flux. Our approach considers this pos- 
sibility by using the best currently available experimental up- 
per limit on the extragalactic proton contribution. We have 
also pointed out in Sect. II C, referring to recent Monte-Carlo 
simulations [15], that the fundamental properties of photohad- 
ronic interactions can affect the bound by up to a factor of 5. 

At neutrino energies below 10 7 GeV, the cosmic ray bound 
computed in this work rises and equals the bound inferred 
from the EGRB flux at about 10 5 GcV, below which the 
EGRB constraint is tighter. In the same energy region, we 
have also shown that effects from extragalactic magnetic fields 
come into play; they could either increase or reduce the bound, 
depending on the details of the field strength and structure, 
and the distribution of cosmic ray/neutrino sources. We there- 
fore do not regard the bound based on the cosmic ray flux as 
a firm upper limit on neutrino fluxes below about 10 5 GeV, 
the main energy region for underwater/ice Cherenkov exper- 
iments such as AMANDA [75]. Models which predict neu- 
trino emission mainly in this energy range are therefore not 
rigorously bounded by cosmic ray data, as for example the 
model of Berezinsky et al. [56] predicting neutrino fluxes 
from cosmic rays stored in clusters of galaxies. Such sources 
could in principle produce neutrino fluxes almost up to the 
EGRB limit within the AMANDA range. However, the neu- 
trino contribution from clusters of galaxies is expected to be 
much lower, as it is limited by their expected contribution to 
the EGRB of -1% [76]. 

At neutrino energies above 10 9 GeV, the upper limit rises 
up to the point where the energy flux of secondary gamma 
rays increases above the level of the observed EGRB. The rea- 
son for the increasing bound is that while cosmic rays from 
evolving extragalactic sources above the nominal GZK cutoff 
reach us exponentially damped due to interactions with the 
microwave background, neutrinos reach us essentially unat- 
tenuated. An observed excess of cosmic ray events above the 
GZK cutoff would therefore correspond to an even more pro- 
nounced excess of neutrinos in the framework of such mod- 
els. Cosmic rays from a fiducial class of extragalactic photo- 
hadronic sources with a sufficiently flat spectrum could thus 
saturate the rising bound at the highest energies and would still 
remain below the local cosmic ray flux at all energies. Clearly, 
this is not the most likely scenario to explain the events be- 
yond the GZK cutoff. For example, if these events are due to a 
single strong nearby source, then we would not expect the ex- 
tragalactic neutrino flux above 10 9 GeV to be at a level higher 
than given by the bound computed by Waxman & Bahcall, 
even if this source would have the spectral properties which 
we used to construct our bound in Sect. Ill C. On the other 
hand, we note that the flux from non-photohadronic sources, 
as for example from the decay of topological defects, can ex- 
ceed this bound because of their different branching ratios 
for the production of baryons relative to mesons [12]. Mea- 
suring the neutrino flux in this energy region, as is planned 
using large air-shower experiments such as the Pierre Auger 
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Observatory [77], would therefore be highly relevant for un- 
derstanding the nature and cosmic distribution of UHECR 
sources. 

The bound we derive for extragalactic neutrinos is indeed 
an upper limit, since we do not consider the possibility that 
ultrarelativistic protons are ejected from their sources without 
preceding isospin flip interactions ("prompt protons")- Obvi- 
ously, any additional prompt proton emission would have the 
effect to lower the relative neutrino flux consistent with the 
upper bound. However, for the most interesting sources, i.e, 
GRBs and AGN jets, the ejection of prompt protons can be 
expected to be strongly suppressed due to the rapid expan- 
sion of the emission region and the implied large adiabatic 
losses. This will allow in the future to use measured neutrino 
fluxes, or experimental limits thereof, as a limit to the con- 
tribution of these sources to the cosmic ray flux. If, for ex- 
ample, GRBs produce the observed UHECR flux and if the 
global GRB emissivity follows a similar evolution as star for- 
mation and AGNs, bursts would have to emit a neutrino flux 
on the level of their optically thin bound, which is roughly two 
orders of magnitude above the original prediction [18]. This 
is independent of the total fractional burst energy converted 
by the protons in their interactions. A factor of ^3— 6 is due 
to evolution, a factor of ~10 due to the limited efficiency of 
p — > n conversion [78], and an additional factor ^5 due to an 
accurate treatment of photoproduction yields in the hard pho- 
ton fields typical for GRBs (note that the yield factors used in 
our figures correspond to AGN-like target photon spectra). In 
this case, the gamma rays would also account for most of the 
EGRB, and one could observe strong TeV bursts from nearby 
sources [79,80]. Currently, only AGN jets are known to con- 
tribute to the EGRB in a major way and this has motivated the 
original model A of Mannheim [10] (which is exactly on the 
level of the optically thin bound as a consequence of the as- 
sumption that AGN jets produce the observed UHECRs). We 
note that the neutrino flux due to nucleons escaping the AGN 
jets and diffusing through the surrounding AGN host galax- 
ies is difficult to assess and limited by the EGRB only [10], 
although this neutrino flux is the one most relevant in the 
1-100 TcV range. 

Generally, the relevance of a bound on optically thin 
sources should not be overestimated. A large number of ex- 
tragalactic neutrino sources could be opaque to the emission 
of UHE cosmic rays producing a neutrino flux well above this 
bound. For the particularly important class of AGN jets, we 
have therefore developed an upper bound to their total neu- 
trino flux using constraints on their gamma-ray and neutron 
opacity set by present gamma ray data. This maximum con- 
tribution of AGN jets to the extragalactic diffuse neutrino flux 
is up to an order of magnitude higher than our optically thin 
bound in the energy range between 10 7 GeV and 10 9 GeV. 
Additionally, we have discussed that for energies below ~ 
10 8 GeV an influence of magnetic fields in the radio lobes of 
active galaxies or on larger scales cannot be excluded. For ex- 
ample, the model of Halzen & Zas [11] predicts a neutrino flux 
about one order of magnitude above our upper limit for AGN 
at E v <~ 10 8 GeV, which is certainly extreme, but still cannot 
be considered as completely ruled out, because at this energy 



the flux of neutrons after turning into cosmic ray protons may 
be diminished by an additional factor ~10 due to adiabatic 
losses in expanding radio halos. At present, the most general 
bound to EGRET detected AGN is given by their contribu- 
tion to the EGRB. The most extreme class of neutrino sources 
opaque to the emission of UHECRs could be AGN jets which 
have not been detected in gamma rays. Prime candidates are 
the GPS and CSS quasars mentioned in Sect. IV A, which to- 
gether make up about 40% of the bright extragalactic radio 
sources. 

Finally, we wish to clarify that our result is not in conflict 
with, but complementary to, the upper limit previously ob- 
tained by Waxman and Bahcall [1]. Their result applies to 
photohadronic sources with a particular spectral shape which 
are transparent to the emission of UHE cosmic rays, such as 
may be representative of low-luminosity BL Lacertae objects 
(like Mrk501) or GRBs. Our result is more general, and there- 
fore less restrictive, indicating that there may be other classes 
of sources, such as quasars, with a different spectral shape, 
and/or which are opaque to the emission of UHE cosmic rays, 
which can produce a higher neutrino flux than the source 
classes considered by Waxman and Bahcall. We also confirm 
the claim by Waxman and Bahcall that large scale magnetic 
fields are unlikely to have an effect for cosmic ray propagation 
at ~ 10 10 GeV, but we additionally point out that magnetic 
field effects cannot be disregarded at lower cosmic ray ener- 
gies. We show that hadronic processes in AGN which are opti- 
cally thick to the emission of UHE cosmic rays could produce 
the extragalactic gamma ray background according to present 
observational constraints. However, the limiting model for 
EGRET-detected AGN presented in this paper is mainly deter- 
mined by the current instrumental limits of gamma-ray astron- 
omy — future observations in the GeV-to-TeV energy range 
may impose stricter bounds, and therefore also limit the pos- 
sible hadronic contribution to the EGRB. Together with an in- 
dependent estimate on the total contribution of AGN jets to 
the EGRB, this may allow one to constrain the overall cosmic 
ray content of AGN jets in the near future. 
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